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Abstract 

The metal-insulator transition in one dimensional fermionic systems with 
long-range interaction is investigated. We have focused on an excitation spec- 
trum by the exact diagonalization technique in sectors with different mo- 
mentum quantum numbers. At rational fillings, we have demonstrated gap 
opening transitions from the Tomonaga-Luttinger liquid to the Mott insula- 
tor associated with a discrete symmetry breaking by changing the interaction 
strength. Finite interaction range is crucial to have the Mott transition at 
a rational filling away from the half filling. It is consistent with the strong 
coupling picture where the Mott gap exists at any rational fillings with suffi- 
ciently strong interaction. The critical regions as a quantum phase transition 
are also investigated numerically. Non-analytic behavior of the Mott gap is 
the characteristic in the weak coupling. It is of the order of the interaction 
in the strong coupling. It implies that the metal-insulator transition of the 
model is of the infinite order as a quantum phase transition at zero tempera- 
ture. Fractal nature of the ground state phase diagram is also revealed. 
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I. INTRODUCTION 



Effects of an electron-electron interaction in electronic systems have become a focus of 
the condensed matter physics recently. In a three or higher dimensions, it is widely believed 
that Landau's fermi liquid theory is valid and the effect of the interaction is absorbed 
into several phenomenological Landau parameters. The system is metallic even with the 
interaction since the ground state is adiabatically connected to the free fermi sea where the 
excitations are given by an electron-hole gapless excitation across the fermi surface. In one 
dimension, however, the ground state is unstable against a perturbation of the interaction 
and the ground state is given by the so-called Tomonaga-Luttinger (TL) liquid. Although 
the fermi liquid theory is not valid in one dimension, the TL liquid is also metallic and 
the excitations are gapless. This TL liquid has been focused recently and there are huge 
number of studies by several techniques as the bosonization and the conformal field theory. 
In the paper, we are trying to investigate the breakdown of the TL liquid behavior in a 
simple fermionic system. In a system with periodic potential, that is, on a lattice as a model 
hamiltonian, the allowed kinetic energy is restricted by a finite band width. Therefore the 
strong electron-electron interaction may bring an opening of the energy gap in electronic 
systems which is known as a Mott transition. The opening of the energy gap implies a 
metal-insulator transition which has become a focus again recently in connection to new 
materials as metal-oxides, organic materials and the high-Tc superconductors. 

As far as the conductivity of the electronic system is concerned, the Mott transition can 
be understood as a freezing of the charge degree of freedom. There can be spin related 
phenomena in a small energy scale, however, the Mott transition is a phenomenon of the 
order of the electron-electron interaction. Noting this fact, we have focused on a charge 
degree of freedom and chosen a model hamiltonian of spinless fermions on a lattice 

H^-tYl + H.c. + J2 Vii-j\ninj, (\4 > 0) (1) 

where the interaction can be long range. ( We set t = 1 in the following.) When the 
interaction is non-zero only for the nearest-neighbor (NN) sites, it is mapped to the spin 1/2 
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antiferromagnetic XXZ model by the Jordan- Wigner transformation. In this case, Haldane 
investigated the model in detail by using the Bethe Ansatz solution of the XXZ model [|I|. 
Generic ground states of the NN model are the TL liquid without the energy gap except at 
a half filling where the model has a metal-insulator transition at Vi = t. It is identified as 
an antiferromagnetic Ising gap in term of the XXZ spin model. 

In this paper, we investigate the model, when the interaction is long range, with various 
filling factors of the fermion numbers. As shown later, it has a rich structure as a fractal. 



II. STRONG COUPLING 

Let us first consider a strong coupling limit V » t of the model when the filling factor 
p = M/L is rational where M is the number of fermions and L is the number of sites. We 
use a periodic boundary condition in the following. Let us assume the interaction satisfies 
the downward convex condition, that is, (z + j)Vi < jVi^i + iVi+j { I < j, « L ) to 

avoid a formation of charge clustering. (See later.) A possible form of the interaction which 
we use in the paper is 

V, = Vf}{a) (2) 

//(«) = (" > 1) (3) 

(-sm7ri)° 

which reduced to a simple power law decaying interaction Vj = V / when j « L in a 
sufficiently large system. The nearest neighbor interaction can be recovered also by taking 
a — >• cxD and L — > oo. When the interaction is sufficiently large, the ground state of the 
system for the rational filling p = p/q with mutually prime integers p and q was known for 
any p and q. It is given by a one dimensional Wigner crystal with period q 0. Although 
the ground state charge ordering ( crystal structure ) is complicated if p and q are large 
integers, it is explicitly given 0. In Fig.|I], shape of the charge ordering is shown for p = p/q 
with g = 113 and p = 1, ■ ■ ■ , g — 1 as an example. As is expected for the formation of 
the charge ordering, there are some commensurate conditions to stabilize the system. This 
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commensurability condition brings a fractal structure into the system as shown later. The 
chemical potential in the thermodynamic limit is evaluated by a similar consideration applied 
for the long range Ising model [Q. It is written as 

oo 

/ioo(p + 0) = y^e+(A;,p) + F(p) (4) 

k=l 
oo 

fiooip-0) = Vj2e4k,p) + F{p) (5) 



k=l 



([^] + l)/[T]-[J]/[T]+i otherwise 



e4A:) = -^]/S^, + (A + l)/[T] 

p ^p' p '-p' 

where [x] is a Gauss symbol which denotes a least integer which is not larger that x and 
F{p) is an order t contribution mainly from the kinetic energy which is a smooth function 
of p. In Figs. 0, the chemical potentials are evaluated for two different interactions. It is a 
devil's staircase which was first discussed by Bak et al. in a context of the long range Ising 
model [0]. The discontinuity of the chemical potential, A/i, which is a key quantity to judge 
whether the system is metallic ( Ap = ) or not (A/i > ) is evaluated as 

P 

A/ioo(p = -) = Poo(p + 0) - poo(p - 0) = Vg{p) 
1 

9iP = ^ = E + f^-i - VT,)- (6) 

This is generically non-negative for the potential which satisfies the convex downward con- 
dition. ( If this condition is not satisfied, A/i can be negative for some filling factor which 
causes an instability. It is identified as a charge clustering.) It implies that the energy gap 
opens for any rational fillings if the interaction is strong enough. Further, as can be seen 
from Eq. (H), the energy gap Ap{p/q) only depends on q. Its dependence is given by a power 
law A/ioo(p = |) ~ '"°^p^^' ■ From the consideration in the strong coupling limit, to have a 
Mott insulator phase in the system with filling factor p = p/q, finite interaction range over 
q sites is crucial. 
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III. WEAK COUPLING 



On the other hand, when the interaction is weak {V << t), the finite band width due to 
the lattice efi^ect ( periodic potential ) is not so important. In this case, one can approximate 
the system as a continuous model with a long range g/r"" interaction and with the periodic 
boundary condition [Q. When a = 2, the continuum model is a Sutherland model which has 
been studied extensively If a = 2, the weak coupling model can be discussed using an 
information from the exact solutions. However, the strong coupling case, of course, neither 
the intermediate coupling case can not be discussed by the exact solutions. The ground 
state of the Sutherland model is a TL liquid without an energy gap independent of the 
filling where the ground state is given by the Jastraw wavefunction. Noting that there is an 
energy gap in the strong coupling, it suggests that there is a finite value of the interaction 
where the energy gap opens. One possibility is that there is alway non-zero energy gap, that 
is, the critical value where the energy gap opens, Vc{p) = 0. However, even in the nearest 
neighbor model where the energy gap can be most stable, there is a gapless TL phase ( XY 
phase in the XXZ model ). Therefore we can expect that there are always finite regions of 
the gapless TL liquid phases in any filling factors and inevitably the Mott transition. In 
the following, we give numerical results to confirm these consideration and study the critical 
region. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

As we have discussed, the ground state of the spinless fermions with long range interaction 
has two phases for any rational filling when the interaction strength is varied. The one is 
the TL liquid metallic phase and the other is the Mott insulator phase. The transition 
between the two is a typical quantum phase transition at zero temperature. In this section, 
these phase transitions are demonstrated numerically. The main focus of the numerical 
calculations here is to investigate a critical behavior near the gap opening ( the transition 
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point ). Our main strategy is to investigate the system from an insulator side. 

We use the exact diagonahzation technique for systems with a periodic boundary con- 
dition. The Hilbert space is divided into several sectors with different momentum quantum 
numbers and diagonalized within them to obtain the lowest few energies. For small systems, 
the full spectra are also obtained. Due to the small system size available, it is not efficient 
to calculate the chemical potential directly. Instead of it, we calculate an excitation gap Eex 
which can be comparable with Afi. 

In Fig. ^ and Fig. ^, lower parts of energy spectra are shown for systems with p = 1/2 
and p = 1/3, respectively. When the interaction is sufficiently weak, one observes a behavior 
of the gapless TL liquid as shown in Fig. ^(a) and Fig. Q^a). On the other hand, opening 
of the energy gap near k = 2nkp, n = 0, ±1,±2, •■■ ( = 27rm/g mod g, m = 1,2, 
for p = p/q) is clearly shown when the interaction is sufficiently strong ( See Fig. ^(b) and 
Fig. §(b)). As is known from the strong coupling analysis, there is a discrete symmetry 
breaking ( translational symmetry ) in the strong coupling phase. For the p = p/q case, 
this is a Zg symmetry breaking. Correspondingly there is almost q degenerate ( ground ) 
states in a finite system. They have different total momentum 2nkF ( n = 1, ■ ■ ■ , g mod q 
) respectively. At these momentum sectors, the lowest energy state is given by one of the 
q degenerate ground states. Therefore the lowest energy gap, the energy difference between 
the lowest energy state at the momentum and the true ground state of the finite system 
( usually k = ), is related to an energy barrier between the degenerate q ground states 
barrier as shown in Eq. (0). The physical energy gap which we concern is the second lowest 
one as seen in the Figures. In the Fig. |^, the lowest energy gap at A; = 2kF for p = 1/3 is 
plotted as a function of the system size for several values of the interaction strength. For 
= 32, the gap size obeys a exponential law as 

Ee.(A; = 2A;f,L) ~e-^^ (7) 

which is a signature of the ( discrete ) symmetry breaking and c is the order of the symmetry 
breaking potential. The discrete symmetry breaking is confirmed numerically for = 32 in 
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Fig. ^. To confirm the discrete symmetry breaking, we have also calculated a spectral flow, 
that is, the energy as a function of the Aharonov-Bohm flux through the periodical system 
( ring) [|^. In Fig. ^, the spectral flows of three different momentum sectors ( A; = 0,27r/3, 
and 47r/3 ) for the p = 1/3 system is shown where the Z3 symmetry breaking is expected. It 
is clearly shown that the 3 low energy states are entangled with each other. That is, these 
three states are equivalent in the thermodynamic limit which is a signature of the discrete 
symmetry breaking. 

Next let us investigate the critical region. For a nearest neighbor model at half filling, 
the behavior of the gap opening is known from the Bethe Ansatz solution and is given by 
A/i ~ exp(— c/ \JV — which is essentially singular at the critical point [^. We have 
numerically investigate the behavior of the gap opening for several filing factors in our long 
range interaction model. In Fig. 0, |, and ^, the excitation gap is plotted as a function of 
the interaction strength. 

In the massless TL liquid phase and the Mott insulator phase with Zq symmetry breaking, 
Ei^{k = 2-71 /L) — Ef^ik = 0) converges to the true energy gap in the thermodynamic limit. 
Therefore we extrapolate Eo{k = 2ti / L) — Eo{k = 0), the energy gap between the lowest 
energy at k = and k = 2tt/L to the infinite size limit by fitting them as a polynomial of 
the 1/L. 



An example of the fitting is shown in Fig. |TO|for p = 1/3 case with several values of the 



interaction strength. The numerical results shown in Fig. 0, |^, and ^ can be well fitted by 
the following essentially singular form as a function of the interaction strength 

EeAV,p) ~ Co(p)exp(-^^^-^^L_), (V > v^ip)). (8) 

It implies that the transition is a generalized Kosterlitz-Thouless (KT) type. Therefore 
the transition is of infinite order as a quantum phase transition. This kind of the singular 
behavior is expected in the conformal field theory for the off-critical behavior near the 
conformally invariant critical phase. In the SU{v) invariant case, the exponent 7 is evaluated 
by the renormalization group analysis to be 7 = i^/(z/ + 2) [Q. At the critical point of the 



model has an apparent Zq symmetry but may have higher symmetry. Then one of the 
possible guesses for the exponent is 7(p = p/q) = q/ {q + 2). 

Numerically it is difficult to determine Vc{p) and 7(p) with sufficient accuracies due to 
the singular behavior of the energy gap. However one may use Vc{p) ~ t/g{p) as a lower 
bound of the Vc since it is plausible to expect the Wigner crystal at the strong coupling melts 
near t ~ Ap^o- This estimate agrees rather well with numerical results obtained (Figs. 0, 
and In Figs. |Tl|, we have shown possible ground state phase diagram using the above 



estimate. In the figure, the interaction strength is plotted in a nonlinear scale. It shows 
that there is alway insulator phase for any rational filling. Since the n — poo curve is the 
devil's staircase, AyUoo(p) shows a fractal structure. Correspondingly, the ground state phase 
diagram of the spinless fermions with long range interaction has a clear fractal structure [jlO[ . 



V. ACKNOWLEDGMENTS 

The author thanks Y. Morita, M. Kohmoto, M. Takahshi, K. Maki and D. Lidsky for 
helpful discussions. This work was supported in part by Grant-in-Aid from the Ministry of 
Education, Science and Culture of Japan. 



8 



FIGURES 

FIG. 1. Ground state charge ordering of the one dimensional spinless fermions with long range 

interaction in the strong coupling. The black points are the positions of the particles, p = p/q 
with q = 113 and p = 1, - ■ ■ ,q — 1. The horizontal axis is the spacial direction j, j = 1, - ■ ■ ,q — \. 

FIG. 2. Evaluation of the chemical potential in the strong coupling limit. The smooth part 
F{p) is set to be zero for simplicity, (a) a = 2, (b) a = 1.5. 

FIG. 3. Low energy spectra of the spinless fermions with long rang interaction classified by the 
momentum for p = 1/2. {a = 2) The different symbols are for different system sizes; A: L = 10, 
O: L = 14, □: L = 18, Q: L = 22. (a) V = 1. The thin lines are guide for eyes, (b) V = 4. 

FIG. 4. Low energy spectrum of the spinless fermions with long rang interaction classified by 
the momentum for p = 1/3. (a = 2, t = 1) The different symbols are for different system sizes; A: 
L = 9, O: L = 15, □: L = 21, Q: L = 27. (a) V = A. The thin lines are guide for the eyes, (b) 
V = 32. 

FIG. 5. The energy gap at = 2/si? is shown as a function of the system sizes for p = 1/3 in 
log scale, {a = 2) The different lines are for different values of the interaction strength; = 4, 8, 
16, and 32 from above. 

FIG. 6. The spectral flows ( energies as a function of the Aharonov-Bohm flux $ through 
the ring system) of the three lowest energy states with momentum k = 0, 27r/3, and 47r/3 for the 
system with the filling factor p = 1/3 {L = 21). {a = 2 and V = 32) $o is the flux quantum. 
When $ = L$o = 21$0) the system returns to the original state by the (small) gauge invariance. 

FIG. 7. The excitation gap Egx as a function of V fov p = 1/2. {a = 2) The data is taken by 
extrapolated to the infinite size. 

FIG. 8. The excitation gap Eg^ as a function of F for p = 1/3. {a = 2) The data is taken by 
extrapolated to the infinite size. 
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FIG. 9. The excitation gap Eex as a function of F for p = 1/4. {a = 2) The data is taken by 
extrapolated to the infinite size hmit. 

FIG. 10. The excitation gap E^xi^L) as a function of L for several values of F ( p = 1/3). 
{a = 2) The solid lines are polynomial fits of the data in 1/L. 

FIG. 11. Estimated ground state phase diagram of the spinless fermions with long rang 
interaction. The solid lines are the region with non-zero energy gap. (a = 2) The rest is a gapless 
TL liquid. The scale in the vertical axis is non- linear as (a) ^ = 1-15 and (b) ^ = 1.01 ( The 

strong coupling region is expanded). 
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